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ABSTRACT 

We analyze the general relativistic oscillations of thin accretion disks around com- 
pact astrophysical objects interacting with the surrounding medium through non- 
gravitational forces. The interaction with the external medium (a thermal bath) is 
modeled via a friction force, and a random force, respectively. The general equations 
describing the stochastically perturbed disks are derived by considering the perturba- 
tions of trajectories of the test particles in equatorial orbits, assumed to move along 
the geodesic lines. By taking into account the presence of a viscous dissipation and 
of a stochastic force we show that the dynamics of the stochastically perturbed disks 
can be formulated in terms of a general relativistic Langevin equation. The stochastic 
energy transport equation is also obtained. The vertical oscillations of the disks in the 
Schwarzschild and Kerr geometries are considered in detail, and they are analyzed by 
numerically integrating the corresponding Langevin equations. The vertical displace- 
ments, velocities and luminosities of the stochastically perturbed disks are explicitly 
obtained for both the Schwarzschild and the Kerr cases. 
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1 INTRODUCTION 

The simplest physical model describing the random mo- 
tion of a non-relativistic Brownian particle is given by the 
Langevin equation, which can be written as (I Coffey et al.l 
12004 ) 



dp 
dt 



-vp + K + Z(t), 



(1) 



where p is the (non-relativistic) momentum, and the three 
terms of the right-hand side of Eq. |l} correspond to the 
friction, to an external force K, and to a stochastic force £, 
respectively, v is a coefficient called the friction coefficient. 
Due to the friction term, the particle in Brownian motion 
loses energy to the medium, but simultaneously gains en- 
ergy from the random kicks of the thermal bath, modeled 
by the random force. The random force, which we assume 
to be a white noise, satisfies the conditions (£(£)} = and 
(£(£)£ = kd (t — t'), respectively, where k is a constant. 
The separation of the force into frictional and random parts 
is merely a phenomenological simplification - microscopi- 
cally, the two forces have the same origin (collision with 
the external medium constituents). The Langevin equation 
has a large range of app lications in physic s, astrophysics, 
engineering, biology etc. (|Coffev et al.ll2004) . 
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An exact formulation of the physics of Barnett relax- 
ation, based on a realistic kinetic model of the relaxation 
mechanism which includes the alignment of the grain an- 
gular momentum in body coordinates by Barnett dissipa- 
tion, disalignment by thermal fluctuations, and coupling of 
the angular momentum to the gas via gas damping was in- 
troduced in lLazarian fc Robergd i|f997l) . The Fokker-Planck 
equation for the measure of internal alignment was solved 
by using numerical integration of the equivalent Langevin 
equation for Brownian rotation. It was also shown that in a 
steady state, energy is transferred back to the grain material 
at an equal rate by Barnett dissipation. 

A simple physical model to describe the dynamics of 
a massive point-like object, such as a black hole, near 
the center of a dens e stellar system was developed in 
IChatteriee et all l)2002h . The total force on the massive body 
can be separated into two independent parts, one of which is 
the slowly varying influence of the aggregate stellar system, 
and the other being the rapidly fluctuating stochastic force 
due to discrete encounters with individual stars. The motion 
of the black hole is similar to that of a Brownian particle in 
a harmonic potential, and its dynamics can be analyzed by 
using an approach based on the Langevin equation. By nu- 
merically solving the Langevin equation one can obtain the 
average values, time-autocorrelation functions, and proba- 
bility distributions of the black hole's position and velocity. 
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By using this model a lower limit on the mass of the black 
hole Sgr A* in the Galactic center can be derived. 

A Langevin-type treatment of the motion of a charged 
particle in the inter-galactic medium magnetic field, which 
allows to estimate both the average and of the root mean 
square t ime-delay fo r part icles of given energy, was intro- 
duced in IVietri et al] (|2003h . This model was compared with 
a scenario where the particles are accelerated at internal 
shocks. The generalized thermodynamics and the collapse 
of a system of self-gravitating Langevin particles exhibit- 
ing anoma lous diffusion in a space of dimension D was con- 
sidered in [Ch avams fc Sird (|2004T ). The equilibrium states 
correspond to polytropic distributions. The dynamical sta- 
bility of stellar systems, gaseous stars and two-dimensional 
vortices was studied by using a thermodynamical analogy. 

There are many processes in the accretion of mate- 
rial onto an active galactic nuclei that occur on relatively 
short timescales as compared to the cooling time of the 
X-ray emitting gas in and around elliptical galaxies and 
galaxy clusters. Accretion rates onto active galactic nuclei 
are likely to be extremely var iable on sho rt timescales. Us- 
ing a Langevin type equation IPopel (|2007l ) has shown that, 
for a simple feedback system, this can induce variability in 
the active galactic nuclei power output that is of much larger 
amplitude, and persists for longer timescales, than the initial 
fluctuations. An implication of this result is that rich galaxy 
clusters are expected to show the largest and longest-lived 
fluctuations. Stochastic variations in the accretion rate also 
mean that the active galactic nuclei injects energy across a 
wide range of timescales. 

Accretion disks around compact objects have become 
standard models for a number of astrophysical phenom- 
ena like X-ray binaries or active galactic nuclei. The disks 
are usually considered as being composed from massive 
test particles that move in the gravitational field of the 
central compact object. Waves and normal-mode oscilla- 
tions in geometrically thin and thick disks around com- 
pact objects have been s tudied exten sively both within 
Newtonian gravity (see iKatd |200ll ) 20 01 for a re- 
view) and within a relativistic framework l|Okazaki et al.l 



_ jw) - 
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disks may produce the quasi-pe r iodic oscillations (QPO's 
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2009) 



in black hole low-mass 
X-ray binaries. 

It is the purpose of the present paper to propose a fully 
general relativistic model for the description of the oscilla- 
tions of the particles in the accretion disks around compact 
objects in contact with an external heat bath. The oscil- 
lations of the disk, assumed to interact with an external 
medium (the heat bath), are described by a general rela- 
tivistic Langevin equation, which is derived, for the case of 
a general axisymmetric gravitational field, by considering 
the perturbations of the geodesic equation describing the 
motion of test particles in stable circular orbits. The energy 
transport in the presence of the stochastic fluctuations of the 
disk is discussed in detail. As an application of our general 
formalism we consider in detail the vertical stochastic oscil- 
lations of the accretion disks in both the static Schwarzschild 



and rotating Kerr geometries. In these cases the equations of 
motion of the disk can be formulated as Langevin equations, 
describing the interaction between the disk and an external 
thermal bath. By numerically integrating the corresponding 
Langevin equation we obtain the vertical displacements, ve- 
locities, and the luminosities of the stochastically oscillating 
disks. 

The present paper is organized as follows. In Section [2] 
we briefly review the motion of the test particles in ax- 
isymmetric gravitational fields. In Section [3] we consider the 
stochastic perturbations of the geodesic trajectories of test 
particles moving around a central compact object. The gen- 
eral system of perturbed equations for stochastic disks is 
obtained in Section [4] The energy transfer in stochastically 
perturbed disks is discussed in Section [5] The vertical mo- 
tions of the disks are considered in Section [6] and it is shown 
that they can be described by a Langevin type equation. 
The disk oscillation frequencies are also obtained for both 
the Schwarzschild and the Kerr metrics. The vertical os- 
cillations of the stochastically perturbed disks are consid- 
ered in Section [7] and the vertical displacements, velocities 
and luminosities are obtained, by numerically integrating 
the corresponding Langevin equation, for both the static 
Schwarzschild and rotating Kerr geometries. We discuss and 
conclude our results in Section [8] 



2 MOTION OF TEST PARTICLES IN 
STATIONARY AXISYMMETRIC 
GRAVITATIONAL FIELDS 

The metric generated by a rotating axisymmetric compact 
general relativistic object can be generally given in cylindri- 
cal coordinates (ct, p, (f>, z) as 

ds 2 = g tt {cdt) + 2g t <t>cdtd(j) + g^^d^ 2 + g pp (dp 2 + dz 2 ) , (2) 

where all the metric functions are functions of p and z only. 
The equatorial plane is placed at z = 0. For a given timelike 
worldline (s) with four velocity u M = dx^/ds and az- 
imuthal angular velocity = d(j>/dt, the specific azimuthal 
angular momentum L of a particle of mass m, given by 

L = Utj, = u (g t j, + g^tt/c) , (3) 

and the specific energy E, given by 

E = -u t = -it* {g t t + gtq>Q/c) , (4) 

are constants of motion. The four- velocity satisfies the condi- 
tion it M it M = —1. In a stationary axisymmetric gravitational 
field the most important types of worldlines correspond to 
spatially circular orbits, given by p = constant, 2 = constant 
and Q — constant, respectively. In this case the four-velocity 
is given by 

u" = w* (1, 0, Q/c, 0) , (5) 
where 

1 



V-gtt - 2g H (fi/c) - g^ {U/cY 



(E-^Y-m 



For circular orbits the four-acceleration a M = 
— {l/2)g a fj^u c "u^ of the test particles has only a radial com- 
ponent, 

1 a B 
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The radial component of the force a p vanishes for two 
particular values Q± of the angular velocity, given by 



o± -gt<p,p±\fi 

C ( 



2 — : 

t<t>,p 



(8) 



3 STOCHASTIC PERTURBATION OF AN 
EQUATORIAL ORBIT 

In the absence of any external perturbation the particles in 
the disk move along the geodesic lines given by 



d 2 x" 
1s2 



+ r 



dx a dx 13 _ 



(9) 



where are the Christoffel symbols associated to the met- 
ric. If the disk is perturbed, a point particle in the nearby 
position a;'' 1 must also satisfy a generalized geodesic equa- 
tion 



d 2 x 



ds 1 



+ r 



>p 



dx'' 



ds 



dx lf) = j_ 
ds m 



(10) 



where / M is the external, non-gravitational force, acting on 
the particle. The coordinates x'^ are given by a/ M = x^+Sx* 1 , 
where Sx^ is a small quantity. Thus, in the first approxima- 
tion, we obtain for the Christoffel symbols 



?&{* + s*) = Kp(*) + Kp,xS*, (ii) 

where A = cT^ a /9x A . Therefore the equation of motion 
for 5x p becomes 



d 2 8x" 
ds 2 



+ 2F 



dx a dSx 13 
" 3 ~ds ds~ 



dr a dr@ i 



jf jm _ q we reo btain the e quation o f pert urbation 
of the geodesic line introduced in IShirokovl (|l973l ) . In the 
following we introduce the four-velocity of the perturbed 
motion as SV = dSx^/ds. <5V M satisfies the condition 
UpSV^ = 0. We consider that the particles in the disk are 
in contact with an isotropic and homogeneous heat bath. 
The interaction of the particles with the bath is described 
by a friction force, and a random force. With respect to a 
comoving observer frame the heat bath has a non-vanishing 
average four- velocity U^. 

In the non-relativistic case the friction force is given by 
fj r — —mvv 1 , where v is the friction coefficient, and v 1 are 
the components of the non-relativistic velocity. The relativis- 
tic generalization of the friction force re quires the introduc- 
tion o f t he fri ction tensor given by (|Dunkel fc Hanggil 
l2005al lbl, 120091 ) 

v£ = um(S£ + 6V a 6V' i ). (13) 
The friction force can be expressed then as 

r Sr = -<{sv a -u a ) = 

-vm(SZ + SV a SVi")(SV a -U a ). (14) 

Let us now introduce a Gaussian stochastic vector field 
\g;x\, where g is the determinant of the metric ten- 
sor, defined by the following correlators: (£ M [g;a;]) = and 
(^[g\xWWM) = D^[g;x-y], where D»» [g-x;y] is the 
noise kernel tensor. From a physical point of view m£ M [g; x] 



represents the stochastic force generated by the interaction 
of the particles of the disk with the external heat bath. Thus 
the motion of the general relativistic massive test particles 
in the stochastically perturbed disk can be described by the 
following general relativistic Langevin type equation 

+ 2 r>W + r^u^sx" = 
-v («g + svasv* 1 ) (sv a - u a ) + e [<?; x] • (15) 

In the following we will use Eq. (|15|) to analyze the phys- 
ical properties of the randomly fluctuating accretion disks. 



4 STOCHASTIC OSCILLATIONS OF 
ACCRETION DISKS 



We assume that the heat bath is at rest, so that (/'' = 
^j^ 1 ' 2 , 0, 0, oj . Also we assume that, since the perturba- 
tions are small, the variation of the gravitational field along 
the z - direction can be neglected. This implies that along 
the equatorial plane gp V ,z = 0, but gp V , Z z 7^ 0. Then from 
the Langevin equation Eq. (|32[1 we obtain the equations of 
motion of the particles in the stochastically perturbed disk 
as 

d 2 5t . „ /„ t „ t Q\ t d8p 



Ftp + r<i„ — 



ds 



-v (<5< + SVcSV*) (8V a -U a )+^ [g; x] , 



(16) 



d 2 8p 
ds 2 



+ 2 ( r « + 0-f + 



tt,p 



t<t>,P n $4>,P 



C ^nr V c 

-v{8 p a + 8V a 5V)8V a +i p [g;x 



(u) Sp- 



in) 



d 2 84> 



2 r 



rl-u 



n 



. dSp 



-v (St + SV^SV 4 ") SV a + [g; x] . 



(18) 



ds 1 



+ 



+ 2F z t 



n 



+ r 



qxb.z 



(u*) Sz 



-u(S z a +SV a SV z )SV a +e [g;x]. 



(19) 



In the absence of the friction force and of the stochastic 
force generated by the interaction between the disk and the 
heat bath we reob tain the equations for th e free oscillations 
of accretion disks l|Semerak fc Zacekll2000l l. 



5 ENERGY TRANSFER IN 

STOCHASTICALLY OSCILLATING DISKS 

In the foll owing we will us e the Eckart model for dissipative 
processes (|Weinberdll972l ). and we will choose a frame de- 
fined by the family of observers moving with the normalized 
four- velocity v p = n M /n parallel to the oscillating matter 
fluid of the disk, where n 2 = — <7 M „n''n". 

We assume that the stochastically oscillating disk con- 
sists of at least one fluid, whose particle number density 
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flux n M = nw M is conserved, nfl = 0, where ; denotes 
the covariant derivative with respect to the metric, and 
n is the particle number density. Let us now introduce 
a Gaussian stochastic tensor field fi mU! ,[<?; a;] defined by 
the correlators (£ M „[gr; as]) = and {^ a p\g;x\^ muv \g;y}) = 
N a B uu\g'-,x\y], respectively, w here (} means statistical aver- 
age ( Hu fc Verdaguerl I2008T 1. The symmetry and positive 
semi-definite property of the noise kernel guarantees that 
the stochastic field tensor ^ v [g;x], or £ M „ for short, is well- 
defined. We assume that this general relativistic stochastic 
tensor describes the fluctuations of the energy-momentum 
tensor ST' 11 ' of the matter in the disk due to the interaction 
with the thermal bath. Hence ST 11 " can be written as 



5T = pvV + 2q^v 



(20) 



where p is the energy density, g M is the transverse momen- 
tum, and 7r M " is the anisotropic stress tensor, as measured by 
an observer moving with the particle flux. g M and -k^ satisfy 
the relations q^v^ = and n^v^ = 0, respectively. Within 
the framework of this physical interpretation we can define 
the internal energy e through the relation n(e + eg) = p , 
where eo is an arbitrary constant (|Hawking & Ellis I \wm . 
The energy-momentum tensor is covariantly conserved, so 
that 5T.£" = 0. By contracting 5T^ V with the observer's 
four velocity w M we obtain 



For the first term in Eq. (|2ip we obtain 



v^ST 



-pv 



(21) 



(22) 



By using the definition of the internal energy and the 
conservation of the particle number density flux we obtain 



Mr"").* = -« - C + v^e v + vptg, 



(23) 



where we have denoted e = e-^v" . The right-hand side in 
Eq. (|21[) can be written as 



(24) 



Thus, the local form of the energy balance of the 
stochastically oscillating disk is given, from the fluid ob- 
server's point of view, by 



(25) 



When the stochastic term in the energy-momentum ten- 
sor is ignored, and by assuming that in this case ST^ takes 
the form ST^u = (p + p) v^v v + pg^ y ', we can define the hy- 
drostatic pressure p of the fluid as the trace of the stress 
tensor 7r M ", p = (l/3)7r?. This allows us to define the vis- 
cous stress tensor P M " as l|Weinberdll972h 



-7T M " + ph! 1 



(26) 



n (e + pit) + qfa + g M i M - v^P^ = 



(28) 



As for the transverse momentum g M it can b e obtained 
from the relativistic analogue of Fourier's law ijWeinberd 
Il972h . namely 



(29) 



where k is the thermal conductivity coefficient and T is the 
equilibrium temperature of the system. Eq. (|28|) is analo- 
gous to the non-relativistic energy balance equation in the 
presence of a stochastic force £, given by 

n^- + V-$- (n-V) •£ (30) 

where e is the internal energy, q is the heat flux, II is 
the total stress tensor, and v is the fluid's three velocity. 
In Eq. (128 p only the term containing the four-acceleration 
lift does not have a Newtonian counterpart. When £ = 0, 
Eq. (|30[) reduces to the standard energy balance equation in 
non-relativistic fluid mechanics ()Landau fc Lifshit j|l987t ). 



6 EQUATION OF MOTION OF VERTICALLY 
PERTURBED ACCRETION DISKS 

In the following we will consider only the vertical oscillations 
of the disk, which, for an arbitrary axisymmetric metric, are 
described by the equation 



at* 2 



+ 



+ 2T$ 



n 



4>,z~ 



+ r 



(pip, z 



(u*) Sz 



-v{8 z a + 5V a 5V z )8V a +r \g;z] 



(31) 



Since <$V" is a small quantity, the equation of motion 
can be further simplified to 

d 2 5z 



l^ + ^ + ul ^ z = i [g ' z] 

where 



(32) 



(33) 



The angular frequency with respect to the radial infinity is 
obtained as 



(?)' 



(34) 



In the following we consider the vertical oscillations of 
the disks around black holes described by the Schwarzschild 
and the Kerr geometry, respectively. 



where is the orthogonal projection operator to the ob- 
server's four- velocity defined as W u = g^ v + v^v v . The vis- 
cous stress tensor can be obtained in a general form as 

P^ = A [a"" + 2^v v) ] , (27) 

where A is the viscosity coefficient and (T M „ = v (y.\v) — 
(l/3)w.™ hpv. By defining the invariant specific volume v as 
the inverse of the particle number density n it follows that 
Wv^ v = vfi = nv. Then for the equation of the local en- 
ergy balance we obtain and 



6.1 Disk oscillation frequencies in the 
Schwarzschild geometry 

The gen eral form of a sta tic axisymmetric metric can be 
given as (|Binni et alj|2005h 

ds 2 = -e 2U c 2 dt 2 + e~ 2U [e 21 (dp 2 + dz 2 ) + p 2 d<t> 2 \ . (35) 
For the Schwarzschild solution 



" 1 - VP^(M + z f + ^ + { M^) 2 -2M 
2 VP 2 + (M + zf + Vp 2 + (M - zf + 2M 
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7 = -In 



VP + zf + p 2 + y/(M - zf + p 2 



AM 2 



4^(M + z) 2 + p 2 ^(M-z) 2 



+ p 2 



-,(37) 



respectively, where M is the mass of the cent ral compact 
object in geometrical units dBinni et al.l 120051 ) . The usual 
Schwarzschild line element in (ct, r, 8, (j>) coordinates is re- 
covered by making the transformations p — \/r 2 — 2Mr sin 8 
and z = (r — M) cos 8. 

For an accretion disk in the Schwarzschild geometry the 
radial frequency at infinity and the proper frequency are 
given by 

4£/-2 7 

Oi (p, z, M) = --U, zz , (38) 

1 - pU, p 

and 



uj 2 ± (p,z,M) = 



•u., 



(39) 



1 - 2pU, p 

respectively. Estimating f2j_ in the equatorial plane at z — 



gives niU =0 = Ml (M + ^M 2 + p< 
obtain 

9 M 



while for u) \ we 



(40) 



p 2 ^jM 2 + p 2 - 2M 2 \ M + yjM 2 + p- 

6.2 Disk oscillation frequencies in the Kerr 
geometry 

The line element of a stationary axisy mmetric spacetime is 
given by the Lewis-Papapetrou metric |Binni et al.| [2005) 

ds 2 = e 2 ^-^ (dp 2 + dz 2 )+p 2 e- 2 ^d4> 2 -e 2i ' {cdt - ujd^) 2 , (41) 

where 7 and ip are functions of z and p only. The metric 
functions generating the Kerr black hole solution in axisym- 
metric form are given by 

1 (R+ + R-) 2 -4M 2 + (a 2 /a 2 ) (R+ - R-) 2 



2 (R++R-+2M) 2 + (a 2 /a 2 )(R + -R„Y 



y = — In 

' 2 



1 {R+ + R-) 2 -4M 2 + (a 2 /a 2 )(R^ 



R- 



m+R~ 

M (R+ +R-+ 2M) [(R+ - R- ) 2 - 4a 2 ] 



; (43) 



(44) 



cr 2 (R+ + R-) 2 -4M 2 + (a 2 /a 2 )(R+- R-) 2 ' 

where M and a are the mass and the specific angular mo- 
mentum of the compact object, R± = \/ p 2 + (z ± <r) 2 and 
a = v M 2 — a 2 . The usual form of the Kerr line ele- 
ment in Boyer-Lindquist coordinates (t,r,6,(j>) is recov- 
ered by performing the coordinates transformation p — 
s/r 2 — 2Mr + a 2 sin 8 and z = (r — M) cos 8, respectively. 

The frequency of the disk oscillations in the Kerr geom- 
etry is given by 

M (zuia + 1373) 



uj ± (p, M, a) 

where we have denoted 

uj! = 4M z/2 R a 5 3 _5 13 12 x 

[8M(2Af + R a )S+ + 2(5M + 3i? CT )p 2 



■0J2 = -SM 3/2 R a SiS™ /z (4MS+ + 2R <y 8+ + p z 



,3 ,19/2 



(45) 



(46) 



(47) 



zu :i = 4R a 8 3 _8% [2(8M 3 - 3MR 2 a - R.l)8 2 + ] + 

4R a 8 3 _8% [3(6M 2 + ZMR a + Rl)8 + p 2 + 3Mp 4 ] , 

(48) 

and 



ZU4 



-4i? CT ^5i° x 



[(8M 3 - 3Mi? 2 - Rl)5+ + 3M(2M + R a )p 2 ] 



(49) 

respectively, with 8± — p 2 + <r 2 ±M and R a = p 2 + a 2 , 
respectively. In the limit a = 0, R a = y M 2 + p 2 and S± = 
R a ± M, and after some algebraic manipulation, we recover 
Eq. (|40[1 giving the oscillation frequency of the Schwarzschild 
disk. 



7 VERTICAL STOCHASTIC OSCILLATIONS 
OF PERTURBED ACCRETION DISKS 

By assuming that the velocity of the perturbed disk oscil- 
lations is small, we can approximate ds as ds m cdt, where 
t is the time coordinate. We also assume that the external 
perturbative force is a function of time only, and neglect 
its possible spatial coordinate dependence. In this case the 
equation of motion of the disk in contact with a heat bath 
follows from Eq. (|32[1 and is given by 



^ + cu^ + c^i (p, M, a) Sz = c 2 ^ (t) , (50) 

where the oscillation frequency of the disk is given by 
Eq. (|40|) in the Schwarzschild case, and by Eq. (|45|) in Kerr 
case, respectively. The random force £ z (t) has a zero mean 
and a variance (£*(*)£*(*')) = (D/c 2 ) S (t - t'), where 
8 (t — t') is the Dirac delta function and D ^ 0. By introduc- 
ing the velocity 8V Z = d8z/8t of the vertical disk oscillation, 
Eq. (|50|) can be written in the form 

d5z = 5V z dt, (51) 
d5V z = -cv5V z dt - c 2 o;i (p, M) Szdt + D 1/2 dW(t), (52) 

where W(t) is a collection of standard Wiener processes. 

A Wiener process W(t) ^ 0, is a one-parameter family 
of Gaussian random variables, with expectations zero and 
covariances E (W(s)W(t)) — min{s,t}. Because the W(t) 
are all Gaussian, this information suffices to determine joint 
probabilities. Alternatively, W(t) may be viewed as a 'ran- 
dom'continuous function with W(0) = 0. A Wiener process 
may be generated at consecutive grid points t n by W(0) = 0, 
W (t n ) = W (t„-i) + (t n — tn-i) Zn, and {Z n } is a set of in- 
dependent standar d Gaussian random var iables (with mean 
and variance 1) (|Skeel fc Izaguirrell2002l ). 

In order to numerically integrate Eqs. (|51[1 an d ([52} 
we cons ider the Briinger - Brooks - Karplus (BBK) 
scheme l|Briinger et al.l [l984). If the time step is taken 
as At, we denote 8z n w Sz(nAt), and we define /„ = 
— c 2 ui\ (p, M, a) Sz„At 2 . The recurrence relation for the 
BBK integrator is 

2 t- l-g/2; 



&z n ^ 



1 + 5/2 
1 

1+772 



1 + 5/2 



Sz n 



+ 



1 + 5/2 



/« + 



DAt 3 ) 



3\l/2 



(53) 
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where we have denoted g = as At. The starting procedure of 
the BBK integrator is given by 

Szi = 8z + (l - f ) V At + i/ + (^^) Z , (54) 

where Szo and Vo are the initial displacement and velocity, 
respectively. The velocity 8V Z can be obtained as 

SV: = SZn ~£ Zn ~ 1 . (55) 
By defining 

E =Kw) 2+ r 2ujl[p ' M > a)dz2 ' (56) 

as the total energy per unit mass of the oscillating disk, 
the luminosity per unit mass L of the disk, representing the 
energy lost by the disk due to viscous dissipation and to the 
presence of the random force, is given by 

(f) a -' 2 ff(«- on 



L = 



dE 
~dt 



7.1 The Schwarzschild stochastic disk 

In the case of an accretion disk in the static Schwarzschild 
geometry of a black hole, the proper frequency of the disk 
oscillations lo\ (M, p) is given by Eq. (|40[) . By expressing the 
radial distance p in terms of the gravitational radius M so 
that p — nM, n — constant ^ 6, we obtain 

1 • (58) 



2 

U}± = 



M 2 [{n 2 - 2) Vl + n 2 



By introducing a set of dimensionless variables (0,5Z), 
defined as t — (M/c) 8 and Sz — MSZ, respectively, and by 
denoting the product vM as f, the equation of motion of 
the stochastic oscillating disk becomes 



d 2 8Z dSZ 
d9 2 +< '~df 



+ 



-5Z = rf (9) . 



(n 2 - 2) s/l + n 2 

where if (9) = M£ z (t). 

The dimensionless luminosity Lg can be written as 



(59) 



dSZ „ 



(60) 



where L e = ML/c 3 . 

The variations of the dimensionless coordinate SZ, of 
the stochastic velocity of the disk, and the luminosity of the 
stochastically perturbed disk are represented, for a central 
object with mass M — lO lo Af0, for n = 7 and n — > oo, and 
for several values of Q, in Figs. []]- [3] respectively. 

In order to describe the stochastic characteristics of the 
physical parameters of the oscillating disk from a quantita- 
tive point of view we compute the Power Spe ctral Distri- 
bution (PSD) of the luminosity ifVa ughan 201ol) (for details 
see the Appendix). The numerical values of the slopes of 
the PSD give an insight into the nature of the mechanism 
leading to the observed variability. The computations of the 
PSD have been done by using the .R software and the as- 
sumption H : P(f) = Pf~ a , where a and ft are constants, 
which was applied to the time series obtained by using the 
BBK integrator. The results for the Schwarzschild case are 
shown in Fig. [4] In the case considered in the present paper, 
the luminosity is determined by the evolution, according to 




2000 4000 6000 8000 10 000 12000 




Figure 1. Variation of the dimensionless coordinate &Z of the 
stochastically oscillating Schwarzschild disk as a function of 8 for 
a central object with mass M = M 10 Mq and Q = 100 (red curve), 
£ = 2.5 X 100 (blue curve), and ( = 5x 100 (green curve). In the 
upper Figure n = 7, while in the lower Figure n — > oo. 



a Brownian-type equation of motion, of a volume element of 
the matter in the disk. For such a situation we expect a PSD 
slope of —2, and this result in fact represents a consistency 
check of the stochastic simulations. 



7.2 The Kerr stochastic disk 



By denoting p = nM, a — kM and 8 = M\Jn 2 + 1 — k 2 , 
with n > and k G [0, 1], the oscillation frequency of the 
Kerr disk, given by Eq. (I45[) . can be represented as 



2 _ cjx (n, fc) 
UJ± ~ M 2 lu 2 (n, fc) ' 

where 
w? (n, fc) 



(61) 



3n +n VTT5 X 
[2k (5 + 38) + 3Vl + 8 (6 + 38 + S 2 )] 



2{l + 8) i/2 x 



4fc(2 + 8) + 



VTTS (-8 + 3S 2 + <5 3 ) 



(62) 



and 



uj 2 2 {n,k) = {1 + 8) [n 2 - 8 (1 + 8) {3 + S)] 2 x |<5- 

[fc + (1 + S) 3/2 ] 2 [4 (1 + 8) + n 2 (3 + 8)] 
[n 2 -8(1 + 8) (3 + 8)] 2 



+ 
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Figure 2. Variation of the dimensionless velocity Vd = d8Z/88 
of the stochastically oscillating Schwarzschild disk as a function 
of 9 for a central object with mass M = M 10 M Q and for ( = 100 
(red curve), ( = 2.5 X 100 (blue curve), and ( = 5 X 100 (green 
curve), respectively. In the upper Figure n = 7, while in the lower 
Figure n — > oo. 



4fc [fc + (l + <5) 3/2 ] 
n 2 -S(l + S) {3 + 8) 



- 1 



(63) 



respectively. 

Then the vertical dimensionless equation of motion of 
the stochastically oscillating Kerr disk is given by 



d 2 SZ JZ 

+ C-T77 + 



dff 2 



dd 



1 (n, k) 
'■ (n, k) 



8Z = rf (0) . 



(64) 



In Figs. [S][7] we have represented the variations of the 
dimensionless coordinate 8Z, of the velocity Vd = d8Z/d9, 
and of the luminosity for a stochastically oscillating Kerr 
disk for a rotating massive central object with mass M = 
1O 1O M0 and a = 0.9M, for n = 7, and for different values 
ofC 

In Fig. [S] we show the PSD of the luminosity of the 
stochastically oscillating disk around a Kerr black hole with 
mass M = 10 10 M Q and a = 0.9M, respectively. 



8 DISCUSSIONS AND FINAL REMARKS 

In the present paper we have introduced a mathematical 
model for the description of the perturbations of a thin ac- 
cretion disk in contact with an exterior stochastic medium 
(a thermal bath). As a result of the interaction between the 
medium and the disk, the motion of the particles in the disk 
has essentially a stochastic nature. To obtain the equations 



Figure 3. Dimensionless luminosity Lg of the stochastically 
oscillating Schwarzschild disk for a central object with mass 
M = M 10 M Q and ( = 100 (red curve), ( = 2.5 X 100 (blue 
curve), and ( = 5 X 100 (green curve), respectively. In the upper 
Figure n = 7, while in the lower Figure n — ¥ oo. 




Frequency 

Figure 4. PSD of the luminosity of the stochastically oscillating 
Schwarzschild disk for a central object of mass M = M 10 Mq. 
In the upper Figure n = 7 and a = 2.076. In the lower Figure 
n — ¥ oo and a = 2.070. The red lines represent the fits according 
to H. 
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Frequency 



Figure 8. PSD of the luminosity for a stochastically oscillating 
disk around a Kerr black hole with mass M = M 10 Mq and a = 
0.9M, for n = 7. The slope of the PSD is 2.046. 



Figure 5. Variation of the dimensionless coordinate &Z of the 
stochastically oscillating disk as a function of 9 for a rotating Kerr 
black hole with mass M = M 1O M and a = 0.9M, for n = 7, 
and for different values of £: £ = 100 (red curve), £ = 2.5 X 100 
(blue curve), and £ = 5 X 100 (green curve), respectively. 




Figure 6. Variation of the dimensionless velocity Vd = dSZ/dd of 
the stochastically oscillating disk as a function of 9 for a rotating 
Kerr black hole with mass M = M 10 M Q and a = 0.9M, for n = 7, 
and for different values of f: f = 100 (red curve), £ = 2.5 X 100 
(blue curve), and f = 5 X 100 (green curve), respectively. 




Figure 7. Variation of the luminosity Lg of the stochastically 
oscillating disk as a function of 8 for a rotating Kerr black hole 
with mass M = M 10 M Q and a = 0.9M, for n = 7, and for 
different values of (: ( = 100 (red curve), ( = 2.5 X 100 (blue 
curve), and ( = 5 X 100 (green curve), respectively. 



of motion of the pert urbed disk we h ave generalized the ap- 
proach introduced in IShirokovl (|l9T3h and lSemerak fc Zacekl 
(2000), by introducing in the equations of the perturbed 
geodesic lines a dissipative viscous term, and a stochastic 
force term, respectively. The equations of the motion of the 
fluctuating disk have been obtained in both the equatorial 
and the vertical plane, as well as the general stochastic en- 
ergy transfer equation. We have studied in detail the vertical 
oscillations of the disk, and we have shown that in this case 
the motion is described by a standard Langevin type differ- 
ential equation for an harmonic oscillator, with a geometry 
dependent oscillation frequency. By numerically integrating 
the equation of motion we have obtained the vertical dis- 
placements, velocities and the luminosity of the disk in the 
cases of both Schwarzschild and Kerr geometries. 

Fluctuations of accretion disks around super-massive 
black holes may have important astrophysical applications. 
The emission of Galactic Black Hole Binaries (BHBs) and 
Active Galactic Nuclei (AGN) displays a significant aperi- 
odic variability on a broad range of time - scales. The Power 
Spectral Density of such variability is generally modeled 
with a power law, P(f) oc f~ a , where power-law index a 
keeps a constant value in a certain range of /, but changes 
among different ranges. At high frequencies, the PSDs of 
both BHBs and AGNs present a steep slope with a ~ 2. On 
the contrary, below a break frequency, typically at a few Hz 
for BHBs, they flatten to a slope with a ~ 1, representing 
the so-called flicker noise (|King et alj|2004 >. However, the 
observed power spectra often deviate from the form / . 
For example, the PSD of Cyg X-l is well described with the 
form Z" 1 in th e soft state, but in the hard state it exhibits 
the form /~ 1,3 jGilfanov Il2010t ). In particular, it was shown 
that the power-law index is around a = 0.8 — 1.3 both in 
the soft state of BHBs a nd in narrow-line Seyfert 1 galaxies 
jjaniuk fc Czernvll2007l ). 

In the case of the present model of the disk oscillations 
under the effect of some stochastic forces the resulting 
Power Spectral Density is of the form P(f) oc /~ 2 . This re- 
sult can be understood qualitatively by noting that P(f) oc 

/ exp (ift) f £ (t') dt' dt oc if' 1 f exp (ift) £(t)dt, which 

gives P oc /~ 2 , since £(t) is a stochastic variable. This re- 
sult is generally also valid for harmonically oscillating, un- 
damped disks, and consequently it is independent of the gen- 
eral relativistic corrections to the equation of motion, since 
these only change the oscillation frequency of the compact 
object - disk system. However, as pointed out previously, 
observations of BHBs and AGNs have shown that the spec- 



© 2002 RAS, MNRAS 000, [L>?? 



Stochastic oscillations of general relativistic disks 9 



tral index is not a — 2, but it is either a w 1, or, more 
generally, a £ (0.8,2). Such a dispersion of the power-law 
index reveals that there must exist some other mechanisms, 
different from purely stochastic oscillations, which may be 
responsible for the observed value of a, like, for example, 
hydrodynamic fluctuations, magnetohydrodynamic turbu- 
lence, magnetic flares, density fluctuations in the corona, or 
variations of the accretion rate, caused by small amplitude 
variations in the viscosity (|King et alj|20()4 ). 

In a model introduced in iTakeuchi fc Mineshigd (|l997l ) 
aperiodic X-ray fluctuations are thought to originate from 
instabilities of the accretion disk around a black hole. To de- 
scribe these ty pe of fluctuations a cel lular automaton model 
was proposed (Mincshig e et a.l.lll994l '). In this model, a gas 
particle is randomly injected into an accretion disk sur- 
rounding a black hole. When the mass density of the disk 
exceeds some critical value at a certain point, an instability 
develops, and the accumulated matter begins to drift inward 
as an avalanche, thereby emitting X-rays. Within the frame- 
work of this model, one can generate / _1 - like fluctuations 
in the X-ray luminosity, in spite of the random mass injec- 
tion. A similar model can be used to explain the optical- 
ultraviolet " flickering" variability in cataclysmic variables 
|Yonehara et al.lll997f ). Light fluctuations are produced by 
occasional flare-like events, and subsequent avalanche flow 
in the accretion disk atmospheres. 

Stochastic oscillations of the accretion disks can also 
provide an alternative model for the explanation of the ob- 
served intra-day variabil ity in BL Lac obje cts, and for other 
similar transient events l|Leung et al.ll20.llh . The basic phys- 
ical idea of this model is that the source of the intra-day vari- 
ability can be related to some stochastic oscillations of the 
disk, triggered by the interaction of the disk with the central 
supermassive black hole, as well as with a background cos- 
mic environment, which perturbs the disk. To explain the 
observed light curve behavior, a model for the stochastic 
oscillations of the disk was developed, by taking into ac- 
count the gravitational interaction with the central object, 
the viscous type damping forces generated in the disk, and a 
stochastic component which describes the interaction with 
the cosmic environment. The stochastic oscillation model 
can reproduce the aperiodic light curves associated with 
transient astronomical phenomena. 

Random or fluctuating phenomena can be found in 
many natural processes. In the present paper we have inves- 
tigated the stochastic properties of the oscillating general 
relativistic thin accretion disks, and we have introduced an 
approach in which the stochastic processes related to disk 
instabilities can be described in a unitary approach by a 
Langevin type equation that includes both the determinis- 
tic, general relativistic, and the random forces acting on the 
accretion disks. Our results can be considered only as a first 
step in the investigation of the stochastic properties of the 
general relativistic accretion disks. In order to obtain a more 
accurate determination of the relevant physical parameters 
of the disks, and to determine more exactly the possible 
astrophysical applications of the model, it would be impor- 
tant to also investigate the radial stochastic oscillations of 
the disks around compact general relativistic objects. 
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series may be expressed in the form of an analytical func- 
tion. For example, we want to test the assumption that 
the source producing the time series behaves so as to be 
described by 

H:P(f) = Pf~ a . 

The software will return the most probable values of the 
parameters {q,/3}, and the probability that the source 
indeed behaves according to this hypothesis. 



APPENDIX 

In statistical analysis, if X is some fluctuating quantity, 
with mean /i and variance a 2 , then a correlation function 
for quantity X is defined as 

= ((X s -y) (X s+T - M )) ^ 
a 2 

where X a is the value of X measured at time s and (} de- 
notes averaging over all values s. Based on the correlation 
function the Power Spectral Distribution (PSD) is defined 



P(f)= / R(r)e 



-i27r/r 



It is straightforward to see the importance of the PSD 
in terms of the "memory" of a given process. For exam- 
ple, if X is the measured luminosity of a perturbed disk, 
the slope of the PSD of a time series of X provides an 
insight into the degree of correlation the underlying phys- 
ical processes have with themselves. The system needs ad- 
ditional energy to fluctuate, and this mechanism is best 
explained for the Brownian motion, in which case the en- 
ergy is thermal. Brownian motion produces a PSD of the 
form P(f) ~ f~ 2 . 

However, assessment of the PSD slope from one single 
observational time-series is non-trivial. In such cases it is 
very helpful to us e dedicated soft ware, such as the statis- 
tical software .R (Vaughan 2010). Basically, the input to 
this software consists of one time series. Any assumption 
about the behaviour of the source that produced this time 
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